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Abstract. — Using the asymmetric fractional calculus of variations, we derive a fractional 
Lagrangian variational formulation of the convection-diffusion equation in the special case of 
constant coefficients. 
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PART I 
INTRODUCTION 

1. Introduction 

The convection-diffusion equation occurs in many physical problems such as porous me- 
dia, engineering, geophysics. It could model the dispersion of a pollutant in a river estuary, 
or groundwater transport, atmospheric pollution, concentration of electron inducing an 
electric current, heat transfer in a heated body. As many PDEs, the solution exists under 
conditions, but is often not known explicitly. Even for linear convection-diffusion equation 
numerical schemes are not always well understood. It is still a challenging problem to ob- 
tain efficient and robust numerical schemes to solve the convection-diffusion equation due 
in particular to the mixing between two different types of behavior, namely the convective 
and diffusive regimes. Let us first recall the equation. Let C M*^ be an open subset with 
a Lipschitz-continuous boundary dQ. The time domain [a,b], < a < b is arbitrary, but 
fixed. Let us consider the general linear parabolic equation of second order: 

ut + il- V)n - div(A' ■ Vn) + l3u = f{t, x) in (a, b] x Q, 

u{t,x) = in{a,b]xdQ, (1) 

u{a, x) = uq{x) in Q , 

where 7 e M*^, A G R'^'"^, /3 G M. 

As an example u is the concentration of a pollutant, transported by a flow of velocity 
7 e M.'^. The tensor A represents the diffusivity of the pollutant specie. The creation or 
destruction of the specie can be taken into account via /3m, and / is the source term. The 
unknown u is both depending on time and space. 

We assume the coefficients are smooth, bounded and satisfying the following properties: 

— A is symmetric, uniformly positive definite s.t. 

A e C°(a, b; L^^inY""^) and 3Ai, A2 > : Aij^l^ < ^A^^ < As^l^ 

— the convection 7 is such that 

^ eC\a,b;W^-°°{QY) and div7 = 

— the reaction /3 is non-negative 

(3 eC\a,b;L°°{n)) a.nd 3/3o /? > /3o • 

These assumptions guarantee that the problem ([T]) is well posed for / G L^{a, b; H^^{Q)), 
and every uq G L'^{Q). 

Let us notice that in the special case of absence of convection, when 7 = 0, the stationary 
convection-diffusion equation is simply the Poisson equation. It is well known that the 
Poisson equation derives from a variational principle also called least-action principle. 
This means that the solution u of the Poisson equation is a minimizer of the following 
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This is not the case of the convection-diffusion equation. For instance, the stationary 
convection-diffusion equation admits the following weak formulation: 



{K-Vv)-V^dx+ {-fV)ucf> + f3u^dx= f4>dx, for any ^ e H^{n). (2) 



Nevertheless the advective term is not symmetric in u and 0. As a consequence, the 
weak formulation ([2]) does not derive from a potential, |19j . This can also be seen as the 
convection-diffusion equation does not satisfy the so-called Helmholtz conditions, i.e. that 
the Frechet derivative of the Euler-Lagrange expression is not self-adjoint. We refer to 
( [14] . Thm. 5.92, p. 364) for more details. 

Let us note that there were some attempt to construct variational formulation for the 
convection-diffusion equation by Ortiz |15] . where he resorts to a local transformation of 
the solution by use of a "dual" problem. In this paper, we prove that the solutions of 
the convection-diffusion equation correspond to critical points of a fractional Lagrangian 
functional. The idea to use fractional derivatives in order to bypass classical obstruction 
to the existence of a Lagrangian functional was discussed by Riewe ( [16] . [T7] ). However, 
for technical reasons his idea cannot provide a fractional variational formulation of the 
convection-diffusion equation. This difficulty was solved recently by Cresson and Inizan in 
[7] . They introduce the asymmetric fractional calculus of variations, and obtain an explicit 
fractional Lagrangian. 

The quest for such a variational formulation is not only an abstract mathematical prob- 
lem in calculus of variations. Due to the existence of different flow regime which depends 
on a coefficient called the Reynolds number, many problems arise in the construction of 
numerical schemes for the convection-diffusion equation. The main difficulty is to avoid 
non physical spurious oscillations which appear in numerical experiment for high value of 
the Reynolds number. Our idea is to use the previous Lagrangian variational formulation 
in order to build a fractional version of a variational integrators. In the classical case (see 
[10] ) the variational integrator possess good numerical properties. A first attempt to con- 
struct fractional variational integrators is done in [4]. 

The outline of the paper is as follow. First we give some notations. In Part [ITl we explain 
Riewe's approach in the case of dissipative equations, the obstructions to the existence of a 
Lagrangian for the convection-diffusion equation and how it is related to the irreversibility 
of the equation. Part IIIII is devoted to recall basics about fractional calculus and introduce 
the asymmetric fractional principle of variations. Finally in Part IIVI we give a Lagrangian 
associated to the convection-diffusion equation based on the fractional derivatives. 
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2. Notations and assumptions 

2.1. Domain. — Let G N. We consider a smooth rf-dimensional bounded convex 
domain Q with boundary dQ. Let {ei, . . . ,ed) be the canonical basis for M.'^. For any 
X G M"^, we denote by the i-th component of x in the canonical basis of M.'^. Let 
1 < i < d and x G f2. We denote by 6i^x the straight line of M*^ defined by 

Si,x = X + Span(ej) 

and Qi^x = ^ Ci 6i^x- As Q is bounded and convex, Qi^x is a segment. Then, it exists 

0'i,x < ^j,a; SUch that 

^i,x ■— {(^1) • • • ) Xi—i, t, Xi^i, . . . , Xd) I t G [ctj (3) 

2.2. Functional sets. — For two sets A and B, J^{A,B) denotes the vector space of 
functions / : A B. Let a,b E M., a < b. Let m,n E N* and p, g G N. Let U be an open 
subset of M*" or the finite interval [a, 6]. The vector space of functions W — j- M" of class 

is denoted by C'p{U). For / G T{[a, b] xQ,R), we denote by: 

X ^f{t,x) t ^f{t,x) 

Let C^''^{[a,b] x fi) and C^{[a,b] x 1]) the functional spaces defined as follow: 

CP''^{[a,b] X n) := {/ G ^([a,6] x 1],M)| Vt G [a,b], f\ G C^^(l]),Vx G fi, G C"'(K6])}, 

andCP([a,6]x(]) := CP'P([a, 6] x fi) when g = p. Let ^^(fi) := {/ G CP(fi) | / = on 
For p = 0, we introduce the following vector spaces: 

C°([a,6]):= {/GC°([a,6]) |/(a) = 0}, 
C°([a,6]) :={/GC°([a,6])|/(6)=0}, 

and for p > 1, let 

C^([a, b]) := {/ G C^'([a, b]) \ f^^\a) = 0, < < |) - 1}, 

C7^([a, b]) := {/ G C*'([a, b]) \ f^'\b) = 0,0<k<p-l}. 
Co^([a,6]) ■.= CU[a,b])nC^[a,b]). 

The set of absolutely continuous functions over [a, b] is denoted by AC {[a, b]) and ACP~^^{[a, b]) 
is the set defined by 

AC^+\[a,b]) := {/ G C*'([a,6]), G AC{[a,b])}. 

Then C^([a, 6]) C AC^([a, 6]). A natural functional space for the study of classical PDEs 
is 

FP''^{[a,b] xQ) := {/GJ^([a,6] xn,M)|VtG [a,b], G C^p, Vx G H, fx e AC\[a,b])}. 

Let FP{[a,b] x fi) := F'''^'([a, 6] x f2) when q = p. For any integer m > 0, H^^iVt) denotes 
the Sobolev space of order m. 
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2.3. Fields. — In this paper we are interested in fields u depending on time t G [a, b] 
and space x & Q: 

u : [a,b]xn — > M . 

(t,x) I — > u(t,x) 

The notation V'u(t, x) e R"^ is the gradient of x i— )■ u{t,x) and dtu{t,x) G M the partial 
derivative of u according to t and dx^f the derivative of / in the i-th space- variable. The 

d 

divergence of a vector field F = (Fi, . . . , Fd) is divF = dx^Fi. Let v : ^l—^'R,l<i<d 



i=l 



and X G f2. We denote by Vi^^ the function defined by 



7 ^ (4) 

y I — > v{xi, . . . ,Xi^i,y,Xi+i, . . . ,Xd). 

li V : f2 — !■ R, we denote also by v its extension to M"^ such that v{x) = if x G R'^\r2. 
For X, ?/ G M'^, X X y denotes the vector in M.'^ defined by 

X X y := (xi?/i, . . . ,Xd|/d)*. (5) 



PART II 

VARIATIONAL PRINCIPLES AND DISSIPATIVE SYSTEMS 



In this section we discuss the problem of constructing a variational principle for dissi- 
pative systems. We review some past issues to solve this problem. Starting from the fact 
that irreversibility was the main obstruction, we propose to use a doubling phase space 
which takes into account the evolution of the system toward past or future and different 
time derivatives operators for each of these variables. We also discuss Riewe's approach to 
dissipative systems using the fractional calculus and prove that it does not give a satisfying 
solution. This part can be avoided in a first lecture and must be considered as an heuris- 
tic support for the mathematical framework of asymmetric fractional calculus of variation 
developed in the next part. 



1. Obstruction to variational principles 

The classical inverse problem of the calculus of variations is the following: Having a 
given set of differential or partial differential equations is it possible to know if they derive 
from a variational principle, i.e. as a critical point of an explicit Lagrangian functional. 
There exists a huge literature on this subject. The Helmholtz's conditions give necessary 
and sufficient conditions for a given set of equations to be obtained as an Euler-Lagrange 
equation. We refer to |14] for more details. Helmholtz's conditions are algebro- analytic 
and are related to the self adjointness of the differential operator attached to the equations. 
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However, they do not give an idea of the physical origin of the obstruction to the existence 
of a variational principle. 

For dissipative systems a classical result of Bauer, [3], in 1931 states that a linear set of 
differential equations with constant coefficients cannot derive from a variational principle. 
The main obstruction is precisely the dissipation of energy which induces a non reversible 
dynamic in time. Bateman, [2], pointed out that this obstruction is only valid if one under- 
stand that the variational principle does not produce additional equations. In particular, 
Bateman constructs a complementary set of equations which enables him to find a varia- 
tional formulation. The main idea behind Bateman's approach is that a dissipative system 
must be seen as physically incomplete. He does not give any additional comments, but we 
think that it is related to the irreversibility of the underlying dynamics. An extension of 
Bateman's construction has been recently given for nonlinear evolution equations [n]. 

Contrary to dissipative systems, for a classical reversible system, the time evolution 
in the two directions (past and future) of time is well described by a single differential 
equation. In fact, we do not know what is the dynamical evolution if we reverse the arrow 
of time for a dissipative system. This is for example the case of the diffusion equation or 
the convection-diffusion equation. Actually, the diffusion phenomenon is non reversible in 
time. 



2. Dealing with irreversibility: doubling phase space and operators 

Bauer's theorem and Bateman construction give two distinct ways to overcome the 
obstruction to a variational formulation: 

— To avoid Bauer's obstruction to the existence of a variational principle, a natural 
idea is to deal with a different kind of differential calculus. This is precisely what we 
will do in the next part, following a previous attempt of Riewe ( [16] . [TT] ) where the 
Riemann-Liouville and Caputo fractional calculus are used. However, as we will see, 
the attempt of Riewe does not work and an appropriate modification of his formalism 
has to be done. 

— The lost of reversibility is interpret by Bateman as an incomplete description of a 
system. We suggest to say that we have two different dynamical variables {x+,X-) 
which represent the evolution of the systems toward future or past. In experiments we 
have only access to the time evolution x+, i.e. that a physical process can be modelled 
by a curve t h-)- x+{t). The evolution of the same physical process if we reverse the 
arrow of time is not known. But this doubling of phase space is only one part of the 
dynamical modelization of a physical process. In experiment, we need an operator 
acting on the phase variables which controls the dynamical behavior. In classical 
mechanics, this is done using the ordinary differential calculus and in particular the 
classical derivative with respect to time d/dt. In the new setting, one must look for 
two time operators d+ and d^ acting on x+ and x_ which we call left and right time 
operators in the sequel. 
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These simple remarks tell us that if we look for a variational formulation of a general 
set of equations then the functional must at least be defined on the doubled phase space 
X = (x+,a;_) with a differential operator D given by DX = (D"x+, where 
and are right and left Riemann-Liouville or Caputo derivatives. 



3. Riewe's approach to dissipative systems: fractional mechanics 

In 1996-1997, Riewe ( [16] . [T7] ) defined a fractional Lagrangian framework to deal with 
dissipative systems. Riewe's theory follows from a simple observation: "// the Lagrangian 



contains terms proportional to ^"^^ J ' ^^^'^ Euler- Lagrange equation will have a term 
proportional to ^ ^ . Hence a frictional force 7^ should derive from a Lagrangian con- 



taining a term proportional to the fractional derivative ^ ^^1/2 j " ^^ere the notation ^^112 
represents formally an operator satisfying the composition rule -jye ° TTE ~ then 

tLb (JjV (Jjb 

studied fractional Lagrangian functional using the left and right Riemann-Liouville frac- 
tional derivatives which satisfy the composition rule formula. Let L a Lagrangian defined 
by: 

L:MxMxMxM — > M, 

(x, u;) I — > L{x,v^,v_,w) 

and X : [a, 6] — )■ M a trajectory. With this notation the functional studied by Riewe is given 

by 



C{x) = [ V]!\, V]/\, ^)dt. 



He proved that critical points x of this functional correspond to the solutions of the gen- 
eralised fractional Euler-Lagrange equation 

d fdL, , / l/2^^ ^ f ^ ^^^l/2^ 



dt \dw 'J ' \dv+'~ 'J ^ ^ 'J dx 

where ^^/^ .- p^^^, I^^^x, f ). 

Riewe derived such a generalised Euler-Lagrange equation for more general functionals 
depending on left and right Riemann-Liouville derivatives and with arbitrary 
a > (see [IB], equation (45) p. 1894). 

1 /2 

The main drawback of this formulation is that the dependence of L with respect to T>jI_ 
(resp. ) induces a derivation with respect to 1)^^ (resp. P^''^ ) in the equation. As a 

l/"^ 1/2 1/2 1/2 

consequence, we will always obtain mixed terms of the form X)_ o Vj^ x or X 
in the associated Euler-Lagrange equation. For example, if we consider the Lagrangian 

L(x,f+,w_,w) = -raw"^ — U{x) + -71'+, 
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we obtain as a generalised Euler-Lagrange equation 



m— — + 7 P_ o X + U [x) = 0. 



However, in general 



^1/2 ^1/2 / CLX 



so that this theory cannot be used in order to provide a variational principle for the linear 
friction problem. This problem of the mixing between the left and right derivatives in the 
fractional calculus of variations is well known (see for example Agrawal [T]). It is due to 
the integration by parts formula which is given for / and g in CQ{[a, b]) by 



In ( |16] . p. 1897) Riewe considered the limit a b while keeping a < b. He then approx- 
imated by . However, this approximation is not justified in general for a large 
class of functions so that Riewe's derivation of a variational principle for the linear friction 
problem is not valid. 

In [5], Cresson tried to overcome this problem by modifying the underlying set of vari- 
ations in the fractional calculus. The set of variations is made of functions h satisfying 
h = h. A critical point of the fractional functional under this restriction is called a 
weak critical point. In that case, we obtain that solutions of the linear friction problem cor- 
responds to weak critical point of the functional but not an equivalence as in the usual case. 
The main problem is that the set of variations is too small to derive a Dubois-Raymond 
result. We refer in particular to the work of Klimek, |12] . for more details. 

In the next section, we review the main result of [7] allowing us to obtain an equivalence 
between solutions of the linear friction problem and critical points of a fractional functional. 



In this part, we recall the classical definitions of the left and right Riemann-Liouville and 
Caputo derivatives in the one dimensional case. We define the multidimensional fractional 
analogous for partial derivatives, gradient and divergence. In particular, we prove a frac- 
tional Green-Riemann theorem. We give a simplified version of the asymmetric fractional 
calculus of variations introduced in [7] for which we refer for more details and weaker as- 
sumptions on functional spaces. This formalism will be used in the last part to derive the 
convection-diffusion equation from a variational principle. 




PART III 



ASYMMETRIC FRACTIONAL CALCULUS OF VARIATIONS 
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1. Fractional calculus 

1.1. Fractional operators: the one-dimensional case. — For a general overview of 
the fractional calculus and more details we refer to the classical book of Samko, Kilbas and 
Marichev, [T5] . 

1.1.1. Fractional integrals. — Let /3 > 0, and /, g : [a,b] — )■ M". 

Definition 1. — The left and right Riemann-Liouville fractional integrals of f are re- 
spectively defined by 

{I'-fm = ^J\r~tr''f{T)dT, 

fort G [a,b], where T is the Gamma function. 

1.1.2. Fractional derivatives. — Let a > 0. Let p G N such that p — 1 < a < p. 

Definition 2. — Let t G [a, 6], the function f is left (resp. right) Riemann-Liouville 
differentiable of order a at t if the quantity 

Kfit)=(^^,°n~'')fit) and 2^^!/(t)= ((-ir^o2^-°^/(t), 

exist respectively. The value V^f{t) (resp. V^f{t)) is called the left (resp. right) 
Riemann-Liouville fractional derivative of order a of f at t. 

Exchanging the order of composition we obtain the left and right Caputo fractional 
derivatives: 

Definition 3. — Let t G [a,b], the function f is left (resp. right) Caputo differentiable 
of order a at t if the quantity 

exist respectively. The value ^'D'^f{t) (resp. '^'Dtf{t)) is called the left (resp. right) 
fractional derivative of order a of f at t. 

In a general setting the subscript refers to Caputo derivative, whereas no subscript is 
used to specify the Riemann-Liouville fractional derivative. 

1.2. Properties. — Let us mention the following relations linking both Riemann-Liouville 
and Caputo fractional derivatives: 
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Lemma 1. — Let < a < 1 and f G AC^{[a,b]) then and exist almost every- 
where and 

V-f{t) = '^Vlf{t)+^j0^f{a) 

i (1 — a) 

The proof can be found in ( |18j . thm.2.2, p. 39). The following lemma concerns the 
semi-group property and the integration by parts formula for the fractional derivatives of 
both Riemann-Liouville and Caputo types: 

Lemma 2. — 1. Composition rule formula: Let f G AC^{[a,b]) , then 

o'^vT = - vT o^vT = -. 

+ + (it + + 

2. Integration by parts formula: Let < a < 1. Let f G AC^([a,6]) and g G Co([a, 6]) 
then 

b rb 



Vlf){t)g{t)dt = I f{t){'V'lg){t)dt 

J a 

VX f){t)g{t)dt = f f{t){^V'lg){t)dt. 

J a 



We refer to [18] . p. 46 for a proof. 

Lemma 3 (Regularity). — Let a > and p G N. 

1. /// G C\[a,b]), then / G Cl{[a,b]). 

2. If f eCl{[a,b]), then V1f = -V1f, and / G ([a, 6]). 

3. If f eCl^\[a,b]), then^Vlf eCl{[a,b]). 

The previous lemma is stated under more general but implicit assumptions (as / G 
X" (L^)) in |18] . However, as we need explicit conditions on the functional spaces in order 
to develop the fractional calculus of variations, we give a less general result but with explicit 
functional spaces. 

Proof — 1. As ^ G C°([a,6]), we have (^)g C°([a,6]) using (Thm. 3.1 of [H], 

p. 53 with A = 0). 

2. It results from 1. and lemma [1] 

3. Let 1 < A; < p. As / G C^+^([a, 6]), then ^ G C^([a, b]) 

d^ d^ df 

__c-T-,a r _ _^7-l-a ^ _ -7-I-Q f(fc+l) 
df^ dt^ + dt~ + ^ ■ 
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As e C%[a,b]), we have J^"" /(fc+i) e using (Thm. 3.1 of [IS], p.53 

with A = 0). Then ^V^f e C^{[a, b]) and -^^T^l / (a) = for 1 < A; < p. Moreover 
■ Vlf{a) = by 1., so that "V^f e C^([a,6]). 



□ 

A similar result holds for the right derivative. 

1.3. Fractional operators: the multidimensional case. — We generalise the previ- 
ous definitions and properties on Q by introducing the multidimensional fractional opera- 
tors. 

1.3.1. Definitions and notations. — For a function v : Q ^ W, we denote by Vi^x the 
function defined on fij .j., as in ([3]) acting in the z-th component of v. The i-th partial 
fractional derivatives are given by: 

the right fractional Riemann-Liouville and Caputo partial derivatives with respect to the 
i-th component of v. In a same way the left fractional Riemann-Liouville and Caputo 
partial derivatives are given by 

d'^v{x) := V^v.^xix,), ^-vix) ■.= 'V'^v,,x{x,). 

The associated Riemann-Liouville and Caputo fractional gradient of u, denoted by 
and '^V"m, are defined by 

/d^u{t,x)\ rd'^u{t,x) 
Vu{t,x):=l : , '=V°M(t,x):= : 

\d2u{t,x)/ Yd2u{t,x) 

The Riemann-Liouville and Caputo fractional divergence of f : f2 — M"' are defined by 

d d 

dW'vix) = J2 dTvi{x) ^dw^'vix) = J2 ""dTviix). 

i=l i=l 

The analogous definitions of the left fractional gradient and divergence hold adding a bar 
on the previous symbols. 

Lemma 4- — Let 7 = (71, ... , 7^) G M^, u G AC^{Q) and x E Q, then we have 

div^/^ ^ ^^^/\(x)) = 7 ■ Vm(x). 

, then 



Proof. — Since 7 x '^V^/'^ 



u 



\ld^d]l^uj 

div^/^ (7 X ^V^'\{x)) = ^^7. (^y' o'^y') u{x). 



i=l 
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Applying lemma [2] concludes the proof. 



□ 



1.3.2. Fractional Green- Riemann formula. — The following lemma is a consequence of 
the unidimensional integration by parts formula of lemma [1] and leads to the fractional 
version of the Green-Riemann theorem mixing Caputo and Riemann-Liouville derivatives 
given in lemma [6l 

Lemma 5. — Let u E AC^{^1) and v E Cq{Q) then 

u{xYd"v{x)dx = / d'^u{x)v{x)dx. 



Proof. — We first extend u, v, ^d'^v and d^u over by associating the value if x G M'^yri. 
In this case, we have 

u{x)^d'^v{x)dx = 



u{x) '^d"v{x)dx 

u{x) '^d'^v{x)dxi 



Xi = — O0 



dxi . . . dxi-idxi+i . . . dxd- 



As VL is convex, there exists a,;..r and hi^x such that Oi.x < hi,x and 



''i,x SI ^i,x 



u{x) '^d^v{x)dxi 



As Ui^x G ^C^{[0'i,xibi^x]) and Vi^x G C'o ([aj,^, the integration by parts formula from 

lemma [2] leads to: 



Ui^xixiYV" Vi^xixi)dxi 



Ui.x{xi)vi^x{xi)dxi 
d°'u{x)vi^xixi)dxi 

+00 

d°'u{x)v{x)dxi. 



As a consequence, we have 
u{x) '^d'^v{x)dxi = 



Xi= — O0 



d^u{x)v{x)dxi 
d^u{x) v{x)dx. 



dxi . . . dxi-idxi+i . . . dxd, 



□ 



Lemma 6 (Fractional Green-Riemann theorem). — Let Q E be a smooth d- 
dimensional bounded convex domain, u E Cl{VL.,M) and v = {vi, . . . ,Vd) E AC^{QY, then 
we have 



v(x) ■ '^V u(x) dx 



div°(f (x)) u{x) dx. 



(6) 
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Proof. — The proof results from lemma O On the canonical basis the scalar product is 
given by 



Applying lemma we have for any 1 < i < d 



n 



d „ 

V"u(x) dx = 2_, / '^d"u{x) dx. 
i=i 



Vi{x)^d^u{x) dx = I d°'vi{x) u{x) dx. 
Q Jq 

The definition of div"" concludes the proof. □ 



2. Asymmetric fractional calculus of variations 

The fractional Euler- Lagrange equations obtained so far in |16L lU [5|, [8] involve both left 
and right fractional derivatives. This is a main drawback, if ones want to recover PDEs 
with order one derivative as composition of fractional derivatives of order 1/2. In this 
section, we give a simplified version of the asymmetric calculus of variations introduced in 
[7] which provides causal fractional Euler-Lagrange equations. We refer to [7] and for 
more details and in particular for weaker assumptions on the functional spaces. 

2.1. Asymmetric fractional Lagrangian. — Here the fractional derivatives are seen 
as partial fractional derivatives according to t and x. With the notations and Ut from 
section II] |2.2I we denote by 

^Vlu{t,x) := ^Vlu^it), ^V"u{t,x) := '=V"nt(x) and Vu{t,x) := Vut{x). 

For a field U = (n^, u_) we define 

"I)°f/(t, x) := ( u+{t, x), - u.{t, x)) ^V'Uit, x) := ("V°M+(t, x), -'=V^M_(t, x)) . 

For a generalised Lagrangian L(t,x,y,v,w, z), we denote dyL, d^L, and dz the partial 
derivatives of L. 

Definition 4- — Let be a Lagrangian L defined as follows: 

L : [a, 6] X X M X R X M"' X M"' — > M 

{t,x,y,v,w,z) I — > L{t,x,y,v,w,z) 

The asymmetric representation of L is denoted by L and given by: 
L : [a, 6] X X M2 X m2 X ]^2d ^ j^2d — ^ ^ 

{t,x,U,V,W,Z) I — > L{t,x,U,V,W,Z) 

where 

L(t,x, {u+,U-), {v+,v-), (2;+, z_)):= L{t,x, u+ + + w+ + w_, 2+ + z_). 

The asymmetric fractional Lagrangian functional £a is defined by: 
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[)3(i+2^ 



Definition 5. — Let L be a Lagrangian as defined in definition^ of class C^([a, h] 
and L its asymmetric representation. The associated asymmetric fractional Lagrangian 
functional of order a is given by 

: C\[a,b] X — ^ R 

U I — > f f L{t,x,U{t,x),''V''U(t),^V'^U{t,x),VU{t,x)) dxdt. 

J a J U 

For convenience, we denote by the functional defined for any U = n_) G 
C\[a,b] X ny by 

La{U){t, x) := L {t, X, U{t, x), '=r'"t/(t)/V"f/(t, x), VU{t, x)) 

= L{t,x,u+{t,x) +u_{t,x), ''V^u+{t,x) - ^V'^u_{t,x), (7) 
^V"u+(t,x) -^V\_(t,a;),VM+(t,a;) + Vu_(t,a;)), 

for any x & Q and t E [a,b] so that 

C^{U)= [ [ L^{U){t,x)dxdt. 
J a Jn 

Let us note that as f/ G C\[a,b] x fi) we have ^^'^(t), '^V^C/(t,x) and V3;?7(t, x) G 
C°([a, fe] X Using the fact that L is Ci([a, b] x M^^+S) obtain that L„(C^([a, 6] x fi) c 
C°([a, 6] X Q) and conclude that is well defined. 

2.2. Asymmetric calculus of variations. — The next lemma explicits the differential 
of the functional defined on U e {C^{[a,b] x Q)y in the direction {CQ{[a,b] x fi))^ 

Lemma 7. — LetU e {C\[a,b]xn)y . Lei*° := {t,x,U{t,x), ^©"^/(t, x)/V°?7(t, x), Vt/(t, x)) 
VKe assume that 

- Vx G fi, t G ACi([a,6]), 

- Vt G [a, 6], X ^ a^L(^") G ACi(n), 

- Vt G [a, 6], X ^ (9,L(*") G Ci(fi). 

r/ien Ca is {C^{[a,b] x VL)Y-differentiable at U and in any direction H = {h+,h^) G 
(CQ([a,6] X Q))"^, the differential of Ca is given by 

DCo,{U,H) = [ /" P^(9^L(*°)+dir(a^L(^"))-div((9,L(*"))l ■/i+(t)dxcit 

+ [ [ [dyL{i.'') - Vld^L{^") - div°(9^L(*°)) - div(9,L(*°))] ■ h_{t)dxdt. 

Proof. — Using a Taylor expansion of the Lagrangian L, we obtain: 

DC^iU,H) = J J [dyLi^)-{h+ + h.) + d,Li*)-CVlh+-'V^h.) 

+d^L{-k) ■ {^Vh+ - ^V"/i_) + d^L{^) ■ (V/i+ + V/i+)] dxdt. 
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As dvL{i^) e y4C^([a, h\) by assumption, using the integration by parts formula of lemma |2] 
with /i+ (resp. in Cl{\a,h]) leads to 

f j d^L{i.)-{^V1h+~^V'^h^)dxdt= [ [ (V^d^Li^) ■ h+ - V^d^L{^) ■ h_) dx dt. 

Ja Jn Ja Jn 

As dwLijk) G AC^iVL) and /i+ (resp. hJ) is in Co([a, 6] x fi), we can apply the fractional 
Green- Riemann formula from lemma to obtain 

/ [ d^L{i<)-{'V''h+-'Wh^)dxdt= I [ (di^''(a^LW)-/i+-div"(9^LW)-/i_)rfa;rft. 

J a Jn Ja Jq^ ' 

The last part of the formula comes from the usual Green- Riemann theorem. This completes 
the proof. □ 

A consequence of the previous lemma is the following characterisation of the extremals 
of La for asymmetric fractional functional as solutions of two fractional Euler-Lagrange 
equations: 

Theorem 1. — Let U E iC^i[a,b] x 

Let := {t, X, U{t, x), "VU^t, x), '=V"t/(t, x), VU{t, x)). We assume that 
-yxen, t^ 9„L(*") G AC^{[a,b]), 

- Vt G [a,6], X ^ d^Liic") G AC\n), 

- Vt G [a,b], X ^ (9,L(*°) G C\Q). 

Then (Co([a,&] x Q))"^ extremals of correspond to solutions of the following set of frac- 
tional Euler-Lagrange equations: 

a„L(*°)+ D-(9,L(^-))+di^"(9^L(^"))-div(9,L(*°)) = 0, 
a,L(^")- P^(a„L(^"))-div"(9^L(^"))-div(9,L(*")) = 0. 

2.3. Specialisation. — This theorem is not applicable as it is for classical PDEs since 
as it provides a system of PDEs. However, by restricting our attention to extremals over 
C^{[a,b] X Q) X {0} over variations in {0} x Co([a, 6] x we obtain a more interesting 
version: 

Theorem 2. — Let u+ G C^{[a,b] x fi). 

Let •k'^ := (t,x,u+(t,x), '^'D" x), '^Vu+(t, x), Vu+(t, x)) . We assume that 

- Vx G fi, t ^ d^,L{^1) G AC^{[a,b]), 

- Vt G [a,6], X d^Lii^l) G AC\Q), 

- Vt G [a, 6], X ^ d,L{i.l) G C^{n). 

Then (m+,0) is a {0} x Co([a,6] x Q)-extremal of the action Ca if and only if u+ satisfies 

dyLi^l) - VI (9,L(^°))-div"(9^L(^^))-div(9,L«))= 0, 
for any x G fi, t G [a, 6] . 

The proof is a consequence of lemma [71 
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Remark 1. — In the previous part, we have heuristically justified the introduction of a 
doubled phase space by saying that irreversibility induces a natural arrow of time. Our idea 
to focus only on curves in C^([a, 6] x f2) x {0} is precisely to say that we are interested 
in one direction of time (here the future). However, the selection of a transverse set for 
variations to the underlying phase space is not so clear. It means heuristically that the 
future depends mostly on the virtual variations in the past. More work are needed in this 
direction. 



PART IV 

APPLICATION TO THE CONVECTION-DIFFUSION EQUATION 



Let us consider the reaction-convection-diffusion equation defined on [a, 6] x by ([T]): 
d 

Truit, x) + 7 ■ Vuit, x) - div (K ■ Vu(t, x)) + I3u(t, x) = fit, x). (8) 
ot 

with constant coefficients 7 G K G W^^'^ and /3 G M. As we aheady mentioned, this 
equation does not derive from a variational principle in the classical sense. Nevertheless the 
result of the previous section allows us to overcome this difficulty and obtain a variational 
formulation of the convection-diffusion equation by mean of the asymmetric fractional 
Lagrangian. Let us defined the extended Lagrangian L given by 

L : [a, 6] X 1] X M X M X X M*^ — > R 

{t,x,y,v,w,z) I — > f{t,x)y--l3y^ + -v^ + -{'yxw)-w--{K-z)-z. 

The direct application of theorem [2] provides that the solutions of the convection-diffusion 
equation are {0} x Cq ([a, b] x r2)-extremals of the following asymmetric fractional functional 
£1/2 defined for U = {u+,U-) by 

Cii2{U)= [ [ L{t,x,u+{t,x) +u.{t,x),^V]^^u+{t,x) -^V^j^u-{t,x), 
Ja Jn 

"V^/'^u+it, x) - u_{t, x), VM+(t, x) + VM_(t, x)) dx dt. 
Namely, the following result holds: 

Theorem 3. — Let u G F'^iVt x [a, 6]). Then u is a solution of the convection- diffusion 
equation (|H]) if and only if {u,0) is a {0} x CQ([a, 6] x Q) critical point of Ci/2- 

Proof. — Let x G and t G [a,b]. Let t G [a,b]. The partial derivatives of L verify: 
- dyL{t, X, u{t, x), ^P^/^ x(t), "V^/2M(t, x), VM(t, x)) = /(t, x) - Pu{t, x). 



dyL{t,x,u{t,x), ^V]!^ x{t),^V^/^u{t,x),Vu{t,x)) = ^V]!^ u{t 



X 



a^L(t,x,M(t,x), ^Vj x{t),^\/^/^u{t,x),Vu{t,x)) = 7 X ^V^/2y(t,x), 
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- d,L{t, X, u{t, x), '^py^ x{t),''V^/^u{t, x), Vu{t, x)) = -K ■ Vu{t, x). 

As ti e F'^{[a,b] X Q) we have that Ux G AC'^{[a,b]) and as a consequence '^V^'^ u G 
AC^{[a,b]). We have also Ut G C^(fi) so that using lemma El we deduce "V^^'^u G C^{fl). 
Moreover Vm G C^([a,6]) so that the conditions of theorem [2] are fulfilled. From lemma 

121 as G AC'^{[a,b]) we have V^"^ o ^V^'^ u = —u. Moreover, as Ut G C^(r2), lemma ID 
applies and we have 

div^/2 X '=V^/^u(x)) = 7 ■ Vu(x). 
This concludes the proof. □ 



PART V 

CONCLUSION AND PERSPECTIVES 



The previous result is only an example of PDEs for which the fractional asymmetric 
calculus of variation provides a Lagrangian variational formulation when this is not pos- 
sible using the classical calculus of variations. As we previously said in the introduction, 
variational formulations of PDEs are important both from the theoretical and practical 
point of view. In that respect our result is far from being complete. 

We have the following list of open problems and perspectives: 

— One must develop the critical point theory associated to our fractional functionals in 
order to provide results about existence and regularity of solutions for these PDEs. 

— Our paper as well as [7] solve the inverse problem of the fractional calculus of varia- 
tions for some classical or fractional PDEs (classical or fractional diffusion equation, 
fractional wave equation, convection-diffusion . . . ). However, we have no characteri- 
zation of PDEs admitting a fractional variational formulation in our setting. In the 
classical case, the Lie approach to ODEs or PDEs as exposed for example in |14] 
provides a necessary criterion known as Helmholtz's conditions. A natural idea is to 
look for the corresponding theory in our case. 

— There exists suitable numerical algorithms to study classical Lagrangian systems 
called variational integrators which are developed for example in |10j . [13] . The 
basic idea of a variational integrator is to preserve this variational structure at the 
discrete level. A natural extension of our work is then to develop variational integra- 
tors adapted to our fractional Lagrangian functionals. A first step in this direction 
has been done in [3] by introducing the notion of discrete embedding of Lagrangian 
systems. However, this work does not cover continuous fractional Lagrangian systems 
and uses only classical discretization of the Riemann-Liouville or Caputo derivative by 
Griinwald-Leitnikov expansions. However for classical functionals we have extended 
this point of view to finite-elements and finite- volumes methods, [6]. We will discuss 
the case of continuous fractional Lagrangian systems in a forthcoming paper. 
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Of course of all these problems are far from being solved for the moment. However, it 
proves that fractional calculus can be useful in a number of classical problems of Analysis 
and in particular for PDEs where classical methods do not provide efficient tools. 
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